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PREFACE 


The work described in this report was carried out under NASA 
Grant No. NAG 3-743 "Interaction Between Tollmien-Schlichting Waves 
and Free-stream Disturbances in Boundary-layer flows". The program 
was monitored by Dr. M.E. Goldstein, Office of the Chief Scientist, 
NASA Lewis Research Center. The principal investigator was Professor 
Richard J. Bodonyi, Department of Aeronautical and Astronautical 
Engineering, The Ohio State University. The results of this study 
will appear in a paper co-authored with J.W.C. Welch, P.W. Duck and 
M. Tad j far in the Journal of Fluid Mechanics. This report 
constitutes the final technical report for this project. 
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SUMMARY 


A numerical study of the generation of Tollmien-Schlichting waves 
due to the interaction between a small freestream disturbance and a 
small localized variation of the surface geometry has been carried 
out using finite difference methods. The nonlinear steady flow is 
of the viscous-inviscid interactive type while the unsteady disturbed 
flow is assumed to be governed by the Navier-Stokes equations 
linearized about this flow. Numerical solutions illustrate the 
growth or decay of the T-S waves generated by the interaction between 
the freestream disturbance and the surface distortion, depending on 
the value of the scaled Strouhal number. An important result of this 
receptivity problem is the numerical determination of the amplitude 
of the Tollmien-Schlichting waves. 


iv 



1. Introduction 


The steady and unsteady effects of small surface-mounted 
obstacles on the boundary-layer flow over a surface have been of 
concern for many years. These effects include most notably 
separation and instability, often leading to transition to 
turbulence. Two main reasons possible for this transition have been 
suggested: either the surface distortion produces, in effect, a 
locally separated shear flow which is susceptible to inviscid 
instabilities associated with the inflectional velocity profile [see 
Smith & Bodonyi 1982, 1985, 1987 and Bodonyi, Smith & Gajjar 1983 for 
a discussion of these inviscid type instabilities], or there is a 
sensitive interaction between the surface distortion and the basic 
flow, possibly with unsteadiness/turbulence in the free stream, which 
can readily accentuate the viscous-inviscid growth of the Tollmien- 
Schlichting instabilities usually present in boundary layers in any 
case. 

The major steady-flow phenomenon observed, separation, is now 
well understood, at least in two-dimensional flows. It is generally 
of an interactive viscous-inviscid type in which the flows inside and 
outside the boundary layer affect each other significantly within a 
relatively short length scale. The question of the stability of the 
separating flow or other locally distorted steady or unsteady motions 
is always present, however, and this has started to receive increased 
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attention, in part because of modern developments in boundary-layer 
methods . 

In this study our concern is with the possibility of the 
generation of Tollmien-Schlichting (T-S) waves due to the interaction 
of small free-stream disturbances and localized variations in the 
surface geometry. In general terms, the sequence of events that 
begins with the excitation of spatially growing T-S waves in a 
boundary layer by free-stream disturbances is known as the 
receptivity problem [Morkovin 1969]. It has been studied by a number 
of authors over the years and the area has been reviewed by Reshotko 
(1976). More recently, Murdock (1980), Goldstein (1983, 1984, 1985), 
Goldstein, Sockol & Sanz (1983) and Goldstein, Leib & Cowley (1987) 
have theoretically investigated the role that small free-stream 
disturbances play in generating T-S waves in boundary-layer flows in 
a variety of situations. In particular we note that Goldstein (1985) 
studied the effect that small variations in surface geometry have on 
scattering weak unsteady free-stream disturbances into Tollmien- 
Schlichting waves. Using the triple-deck scalings of Stewartson 
(1969) , Goldstein concluded that relatively small surface variations 
which provoke correspondingly small pressure changes can produce a 
large coupling between the T-S waves and the imposed disturbance when 
these variations are sufficiently rapid, i.e., when they occur on the 
scale of a T-S wave. Goldstein's analysis provided a qualitative 
explanation of the Leehey & Shapiro (1979) boundary-layer receptivity 
measurements. Further comparisons with the Leehey & Shapiro 
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experiments using a solution of the local Orr-Sommerfeld equation 
have been made by Goldstein & Hultgren (1987) . 

However, Goldstein's analysis is limited in that he took for the 
steady flow the linearized solution of Stewartson (1970, 1971) for 
the interactive flow in the vicinity of a sharp corner in an 
incompressible flow. Thus Goldstein could not consider the effects 
of surface variations of sufficient size to provoke boundary-layer 
separation or even a nonlinear response in the steady flow, although 
his analysis does account for nonparallel flow effects on the 
stability of the flow. 

Since our interest is in the stability of nonparallel flows and 
especially those with strong local streamwise variations in surface 
geometry, it is appropriate to take the steady nonlinear viscous- 
inviscid interactive solutions, of the triple-deck and similar kinds, 
for the basic steady motion. This is because, as is now well-known, 
flow reversal for small or large scale separations occuring in such 
flows is not necessarily a catastrophic event: the solution at the 
separation point is regular due to the presence of interaction, 
unlike that in steady classical, i.e. noninteractive, boundary layers 
for instance. Hence a steady nonparallel basic flow with a small 
localized region of reversed flow can be described fully by the 
classical boundary-layer equations, subject to an unknown pressure 
which must also be computed as part of the solution. Additionally, 
we shall assume that the unsteady flow is governed by the linearized 
Navier-Stokes equations, as discussed in the following section. 
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2 . Problem Formulation 


We wish to study the interaction between an unsteady freestream 
and a small surface perturbation on a flat plate, such as a hump or 
trough, for an incompressible two-dimensional viscous flow. Thus 
following Goldstein (1985) we take the upstream motion to consist of 
a uniform flow with velocity U* plus a small harmonic perturbation of 
frequency 0 and constant amplitude u* « U*, so that the unsteady 
motion can be analyzed as a linear perturbation of the uniform steady 
flow, U , i.e., U (1 + u /U e ). 

Consider a Cartesian coordinate system (x*,y*) with x* tangent to 
and y* normal to the flat plate with the origin taken at the leading 
edge. Further, define the Reynolds number Re = u*L*/* / # where L* is 
the distance of the surface perturbation from the leading edge of the 
flat plate and v is the kinematic viscosity of the fluid. For 
convenience, we introduce the small parameter c ■ Re - */ 8 and consider 
solutions of the Navier-Stokes equations when Re » 1. Specifically, 
we wish to consider the problem of flow over a small hump of length 
0(L*e 3 ) and height 0(L*c 5 ) with a profile of the form 

y /L - c 5 F({X -L }/e 3 L ) , (2.1) 

with F = 0(1), positioned at a distance L* from the leading edge. The 
interaction region is shown schematically in Figure 1. 

Following Goldstein (1985) we write 
U(x,y,t) = U Q (x,y) + $u(x,y,t) , 

V(x,y,t) - V Q (x,y) + 5v(x,y,t) , 

P(x,y,t) * P Q ( x , y ) + *p(x,y,t) 
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( 2 . 2 ) 

(2.3) 

(2.4) 



are 


where the steady velocity components , U Q , V Q , and pressure, P Q 

* *2 - - - 
normalized by and p , respectively and u,v,p the unsteady 

velocity and pressure terms normalized by u* and pu*U*, 

respectively. Furthermore, we define 

x * (x*-l*)/l*, y = y*/L* t t = nt*, 6 = u*/u*. (2.5) 

GO GO 

Substituting (2.2) - (2.4) into the Navier-Stokes equations and 
neglecting terms of 0(5 2 ), we get the linearized perturbation 
equations 
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(2.8) 

where 




S = QL*/U* 

' CO 

/ 


(2.9) 


is the Strouhal number. 

Finally, we note that the physical interaction between the 
oncoming boundary layer, freestream disturbance, and hump is governed 
by a triple-deck structure, centered near the surface distortion. 

The details of the structure as applied to this problem have been 
given by Goldstein (1985) and we, therefore, only summarize the 
relevant portions here. As is usually the case the viscous 
interaction problem essentially reduces to a study of the lower-deck 
equations. Thus the appropriately scaled variables in the lower deck 


for the steady flow are 

U Q (x,y) - cU(X,Y) + 0(c 2 ) , (2.10) 

V Q (X,y) - c 3 V (X, Y) + 0(c 4 ) , (2.11) 

P 0 (x,y) * c 2 P(X) + 0 (c 3 ) , (2.12) 

where X « x/e 3 , Y * y/c 5 . (2.13) 
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U, V, and P are found from the solution of the lower-deck equations: 


u x + v y - 0, (2.14) 

mJ x + VU y = -P x + U yY , (2.15) 

subject to the boundary conditions 

U = V = 0 on Y = F (X) , (2.16) 

(U,V,P) - (Y,0,0) as |Xj - « , (2.17) 

U - [Y + A(X) ] , Y - • , all X, (2.18) 


where we have assumed that a simple renormalization of the variables 
has been carried out in order to set the value of the wall shear of 
the oncoming, undisturbed boundary layer to unity. 

Finally, the interaction condition for incompressible flow is given 
by the Cauchy Hilbert integral 

P (X) = lAj_"(X-$)“ 1 (dA/d^)dC . (2.19) 

This steady problem was originally formulated by Smith (1973). 

For the unsteady flow we introduce the following lower- 
deck variables 

u(x,y,t) * e” it u(X, Y) + 0(c) , (2.20) 

v(x,y,t) - c 2 e” lt v(X,Y) + 0(c 3 ) , (2.21) 

p(x,y,t) « ce -it p(X) + 0(c 2 ) . (2.22) 

In these expressions we have utilized the fact that since the 
unsteady flow is governed by the linearized Navier-Stokes equations 
we can seek solutions which have a harmonic time dependence. 
Substituting (2.13), (2.20) - (2.22) into (2.6) - (2.8) yields, to 
leading order in c , 

-ic 2 Su + Uu x + uU x + Vu y + vU Y = 

6 


" P X + U YY ' 


(2.23) 


(2.24) 


U X + V Y ~ 0 

As noted in the Introduction, our interest in this paper is in 
the relatively high frequency case where we choose 0 to be of the 
same order as the Tollmien-Schlichting wave frequency at and upstream 
of the lower branch of the neutral stability curve. For this reason 
we require that 0 = 0(c 2 ) i.e., the Strouhal number S = 0 (e“ 2 ) . Thus 


we define a scaled Strouhal number S such that 

o 


c 2 S, 


0 ( 1 ). 


(2.25) 


o ' o 

Furthermore, in this case the Stokes-layer thickness is found to be 
of 0(c) also, and, therefore, it will be of 0(c 5 ) in terms of the 
lower-deck scalings. Hence the Stokes-layer thickness is the same as 
that of the lower deck. Upstream of the triple-deck region, where 
the mean flow changes on the scale of x, the unsteady flow in the 
boundary layer is given by the Stokes solution, which can be written 
in terms of the lower-deck scalings, as 

u = i - , (2 . 26) 

P * iXS o . (2.27) 

The lower-deck problem is completed by solving (2.23), (2.24) 
using (2.25) subject to the no slip condition at the wall 

u-v=0 on Y * F(X) , (2.28) 

matching with the main-deck solution 

u - 1 + a (X) , Y - all X, (2.29) 

and matching with the upstream Stokes layer solution given by (2.26) 
and (2.27) for X -* -®. Finally, the relationship between a(X) and 
dp/dX is given by the interaction condition 


d 2 a/dX 2 = l/*J_“[dp/d£-iS o ]/(*-X)d£ . 


(2.30) 
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It was found that numerical solutions of the disturbance 
equations were not obtainable for supercritical disturbances using 
the interaction law (2.30). Thus an alternative method used by 
Bodonyi & Duck (1985) in treating three-dimensional interacting flows 
has been employed here. In this approach, the relationship between 
the pressure p(X) and displacement thickness -a(X) is found through a 
numerical solution of the upper-deck equations as opposed to the 
Hilbert integral representation (2.30). Specifically, it can be 
shown that the appropriate boundary-value problem in the upper-deck 


for the disturbance pressure is given by 

PXX + Pyy = 0 ' (2.31) 

with boundary conditions 

Pa(x, 0) = d 2 a/dX 2 , (2.32) 

P(X,y) -> 0, y - «, all X , (2.33) 

P(X,y) + iS Q X - p(X) as y -* 0, (2.34) 

p - 0 as as X -* -«, (2.35) 

P x ~ ikp - 0, as X - ®, (2.36) 

where we've written 

P - P. “ « tP(X,y) + is o X] , y - y A 3 . (2.37) 


Note that (2.36) defines a radiation condition applied on the 
disturbance pressure at the downstream boundary to simulate the 
outward propagating pressure disturbances there. The wave number, k, 
which depends on S Q , is found from the solution of the classical Orr- 
Sommerfeld eigenvalue problem for the Stokes' layer flow. 
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Alternatively, k can be computed iteratively from the numerical 
computations, as will be discussed below. 

The entire unsteady solution is thus known once the solution of 
this boundary-value problem has been obtained. As a special case 
Goldstein (1985) gives an analytical solution to the problem using a 
linearized solution of the steady-flow problem for a slightly cranked 
flap on a flat plate due to Stewartson (1970,1971). Utilizing these 
results for U(X,Y) , V(X,Y) and P(X) Goldstein was able to solve the 
unsteady problem using a Fourier transform technique. 

In this paper we solve the corresponding problem for surface 
distortions such that the steady flow is fully nonlinear in 
character. However, in this case both the steady and the unsteady 
problems must be solved numerically. Before discussing the numerical 
approach used in this study it is appropriate to note that the 
problem formulated above differs in a significant way from the usual 
hydrodynamic stability problem which generally involves the solution 
of the Orr-Sommerfeld equation. The Orr-Sommerfeld eigenvalue 
problem when applied to boundary- layer flows leads to the solution 
for "free" disturbances which are, in fact, the normal modes of the 
boundary layer and they are usually referred to as the Tollmien- 
Schlichting waves. While important in their own right, the normal 
mode representation of a small disturbance spectrum cannot be 
conveniently extended to finite-amplitude, i.e., nonlinear 
disturbances, nor are they useful by themselves in understanding how 
external disturbances, such as free-stream turbulence or surface 
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roughness, etc, feed energy into the boundary layer, thereby exciting 
spatially growing T-S waves. 

One approach to understanding both the linear and nonlinear 
evolution of T-S waves in a boundary layer is to consider an initial- 
value problem, using the unsteady triple-deck equations. Such 
studies have been carried out recently by Smith (1985) , Smith & 
Burggraf (1985), and Duck (1985). 

Alternatively, one can study the nature of the coupling between 
an imposed free-stream disturbance and the growth of T-S waves in the 
boundary layer. This receptivity problem differs from the 
hydrodynamic stability problem both physically and mathematically. 
Physically, it is the response in the boundary layer to some 
externally imposed disturbance. Mathematically, it is no longer an 
eigenvalue problem. Instead it is a boundary-value problem as can be 
seen from equations (2.23) - (2.24) and (2.26) - (2.30), wherein the 
boundary layer is driven by some external forced oscillation with its 
response being a solution of the linearized disturbance equations 
having the same frequency and phase speed as the particular forcing 
disturbance being studied. The primary objective of this work is a 
detailed numerical study of this receptivity problem for a range of 
values of S Q , which represents the nature of the freestream 
disturbance, and a representative surface distortion which we take to 
be 

F (X) - h(l + X 2 )" 1 , (2.38) 

where h is an order one factor which gives the height of the 
distortion relative to the lower-deck scalings. 
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3. Numerical Method 

3.1 Steady-State Solution 

First we consider the steady flow problem defined by equations 
(2.14) - (2.19) along with (2.38). Numerical solutions have been 
found using a finite-difference procedure developed by Smith & 

Bodonyi (1985). Briefly, the governing equations are replaced by 
difference representations for U * r = Uy, and P with uniform 
steps in X, Y. The computational domain extends from X * X 1 (<0) to 
X - X 2 (>0) and from Y * 0 to Y = Y #( with starting conditions (2.17) 
specified, in effect, at X = X^ 

The nonlinear difference equations at a given streamwise location 
X are solved to within a tolerance of 10 in absolute value by 
Newtonian iteration using Gaussian elimination and back substitution. 
The solution is then advanced to the next streamwise location and the 
process repeated until the entire domain is covered. Since the 
problem is interactive, multiple forward-marching sweeps are 
necessary until a tolerance of 10 between successive values 
obtained for P(X) is satisfied for all X. At this point the solution 
is said to have converged in the global sense. The diagonally 
dominant nature of the finite-difference form of the interaction law 
(2.19) makes this multi-sweeping process both fast and stable. 
Whenever flow reversal occurs, i.e., U < 0, windward differencing is 
used to represent UU X in finite-difference form. 

The numerical solution has been found using (2.38) for the 
surface shape with h * 0.1, 1.0 and 5.0. Representative 
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distributions of the wall shear r (X, 0) , and pressure, P(X) , of the 
steady flow are given in Figure 2 for h = 5.0. 

With the steady solution known, we now proceed to consider the 
numerical solution of the complex unsteady linearized boundary-layer 
equations (2.23) - (2.24). In our approach, the governing equations 
were solved in the physical plane using a finite difference method. 

3.2 Unsteady Solution 

The unsteady equations (2.23) - (2.24) were replaced, after 
subtracting out the Stokes shear-wave solution, by a system of 
difference equations of second-order accuracy to be consistent with 
the numerical method used for the steady flow problem. Since the 
governing equations are linear no iteration in the normal direction 
is necessary at a fixed streamwise location. A single sweep across 
the boundary-layer region was sufficient to determine the solution 
there. Thus one complete sweep of the computational domain could be 
accomplished quickly. Multiple sweeps of the entire domain are still 
necessary to obtain the global solution, however, due to the elliptic 
nature of the interaction law (2.30) or equations (2.31) - (2.36). 

Initial attempts to solve the problem were made using the 
pressure-displacement interaction law (2.30) in a form utilizing the 
ideas first put forward by Veldman (1979) and fully discussed by 
Smith & Bodonyi (1985) . Indeed, for values of S Q < 1 acceptable 
solutions could be found. However, as the scaled frequency, S Q was 
increased towards its critical (i.e. neutral) value of S 0cr ^ t ® 2.296 
acceptable numerical solutions became increasingly more difficult to 
obtain. It appears that these difficulties are related to the use of 
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the interaction condition (2.30). Numerically, the Hilbert integral 
is truncated to the finite range X 1 < £ < X 2 , thus it is implicitly 
assumed that the tails of the integral over -® < £ < X x and £ > X 2 
are negligibly small. This assumption appears to be satisfactory for 
the upstream "tail", but it is questionable for the downstream 
behavior of the integral. Furthermore, in addition to the 
difficulties associated with the finite integration domain, there is 
some question as to the convergence properties of the Cauchy-Hilbert 
integral itself, at least when the disturbances are supercritical, 
due to the exponentially growing form of the disturbance quantities 
downstream . 

To overcome these difficulties, the method developed by Bodonyi & 
Duck (1988) for solving flows with viscous-inviscid interaction has 
been succesfully applied to this problem. In this method the 
relationship between the disturbance pressure and displacement 
thickness is found through a numerical solution of the appropriate 
upper-deck equations (2.31) - (2.36) as opposed to the Hilbert 
integral representation of the solution, equation (2.30). The 
crucial feature of the scheme is the inherent numerical coupling 
between the viscous boundary-layer solution and the inviscid outer 
flow solution which is carried out simultaneously in the spirit of 
the scheme proposed by Veldman (1979) . Using this approach, the 
difficulties associated with the convergence of the Hilbert integral 
for supercritical disturbances can be avoided and the proper 
downstream boundary condition (radiation condition, see 2.36) can be 
applied to the disturbance pressure. 
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For numerical convenience we apply the Prandtl transposition law 


and we also subtract out the Stokes solution. Thus consider the 


following change of variables 

u(X,Y) = 1 - e 1 ° + u q (X,Y) , 

v(X, Y) * v q (X,Y) , dp/dX - iS Q - dp^dX . 
Then the disturbance equations can be written as 

u oX + v oY " 0 ' 

(D x - is 0 )u 0 + au oX + U Y v o + vu oV + d Po /dx - u oiY 
-U x + tu x ♦ 1 3/2 s l/2 v)e iV2s i / 2Y # 


with boundary conditions 


u_ = v * ° 
o o 


on Y = 0, all X , 


u -* 0 as X -* -® , 
o 


u_ - a(X) + e 
o 


all Y , 

i 3 /2 s £/2y 


as Y all X. 


Given a guess or an update for u Q , v Q , P Q (X) and p(X,y) 


(3.1) 

(3.2) 

(3.3) 

(3.4) 

(3.5) 

(3.6) 

(3.7) 


everywhere, equations (3.3) - (3.7) are inarched forward in X, while 
simultaneously equations (2.31) - (2.36) are solved along a line of 
varying y. This then determines the complex-valued functions u Q , v Q , 
p and p (and hence a(X)) at a given streamwise location X. Sweeping 


through all X stations constitutes one global iteration. Convergence 
is finally attained when a global convergence test on the disturbance 
displacement thickness -a(X) is satisfied. 

The main features of the numerical scheme are the following. Two 
and three-point differencing in Y is used for equations (3.3) and 
(3.4), respectively, wih equation (3.7) applied at Y - Y^. Three- 
point central differencing is used to approximate (2.31) in both 
dimensions, while condition (2.33) is applied at y = y o . Equation 
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(2.32) is approximated by one-sided differencing in y and a second- 
order scheme for X derivatives. Finally, the radiation condition is 
applied in the following form to estimate the disturbance pressure at 
the downstream boundary 

P^max'*) = tP^ X max“ 2AX ^ ) ]/[2ikAX-3] . 

The value for k is either prescribed as discussed earlier, or by 
estimating its value from the relation (3 p/d X) /ip from values of X 
reasonably far downstream, and then feeding this value back into the 
numerical computations. Numerically, the results indicate only 
slight differences in estimating k in these two ways. This isn't too 
surprising since the viscous-inviscid interaction is a local 
phenomenon and the behavior far downstream should approach that of 
the classical stability theory. 

Supposing we have n points in Y and m points in y, then at each X 
station, the difference approximation of equations (3.3), (3.4) and 
(2.31), together with the interface conditions (2.34) and (3.7) can 
be written conveniently in matrix form as discussed by Bodonyi & Duck 
(1988) . The overall scheme is nominally second-order accurate in the 
grid spacings AX, AY and Ay. The resulting matrix equation is then 
solved using standard Gaussian elimination procedures and back 
substitution. 
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4. Numerical Results 


In this section the results of the numerical computations will be 
presented. As part of a collaborative effort to confirm the finite- 
difference numerical solutions. Dr. P.W. Duck and Mr. J.W.C. Welch of 
the University of Manchester, Manchester, England, have solved the 
same receptivity problem discussed in the previous section using both 
time-marching and time-periodic spectral methods. These spectral 
computations were not supported by this NASA grant. However, in the 
figures that follow, the finite-difference results are represented by 
solid curves and the spectral results by the symbol o for comparison 
purposes. A detailed discussion concerning the spectral method as 
applied to this problem can be found in Bodonyi, Welch, Duck and 
Tad j far (1989). 

We first consider the effect of the Strouhal number, S Q , on the 

interaction between the unsteady flow and the surface distortion. To 

minimize any nonlinear effects, a small value of the hump height, h = 

0.1, for which separation does not occur in the mean flow, was 

chosen. The numerical solution was found for several values of S , 

o 

ranging from 0.5 to 3.5. In most cases 200 points were taken in the 
streamwise direction over the range -10 s X s 10 (AX»0.10). To solve 
the lower-deck equations, 50 points were taken across the lower-deck 
region over the range 0 s Y s 7.5. Also, an* additional 50 points 
were used in the upper-deck y scaling over the range 0 s y s 5, to 
solve Laplace's equation for the pressure in the upper-deck region. 

To test the sensitivity of the thickness of the upper-deck region, in 
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some cases the upper-deck range was extended to y = 10 (100 points), 
but no significant differences in the results were noticed. 

Using the finite difference method of solution, convergence of 

the numerical computations was achieved when the absolute value of 

the difference in the displacement thickness, -a(X), between two 

successive iterates was less than 10 ^ for all X. The number of 

iterations required for convergence was found to be quite sensitive 

on the value of S under consideration and also on the initial 

o 

guesses taken for the disturbance profiles. For a subcritical 
disturbance of S Q = 0.5, less than one hundred iterations were 
sufficient for convergence. However, for a supercritical disturbance 
with S Q » 3.0 approximately 3400 iterations were necessary to obtain 
converged results. Also, by gradually increasing the value of S Q and 
using the converged solution for the previous value of S Q as the 
initial guess, we could reduce the number of iterations required and 
also maintain a small maximum fluctuation between successive 
iterative values of the displacement thickness. 

The disturbance produced by the interaction with the hump should 
ultimately decay sufficiently far downstream of the hump if the 
scaled freestream Strouhal number is below its critical value, S 0crit 
» 2.296. To illustrate this eventual damping of the unsteady 
disturbances, an extended range of -10 s X s 40 (aX= 0.25) was 
considered. The real parts of the complex-valued disturbance 
pressure and wall shear distributions for a representative 
subcritical case, S Q » 2, are presented in Figure 3 from both the 
finite difference and spectral methods. The decay in the disturbance 
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amplitude for all quantities is clearly seen for X a 10. Similarily, 
the disturbances should amplify downstream of the surface distortion 
if S Q is supercritical. The same extended X range was thus 
considered for a typical supercritical case, S Q = 2.5 and Figure 4 
shows the amplification of the disturbances in the streamwise 
direction in this case. 

The supercritical results given in Figure 4 were calculated 
using a linearized time-marching scheme and are shown for t = 14*. 

For both cases, excellent agreement exists between the finite 
difference and spectral method results for the positioning of the T-S 
waves. There is a discrepancy of approximately 10% in the wave 
amplitudes. This difference can be reduced by increasing the number 
of grid points in the X direction of the finite difference method as 
discussed below. 

Numerical experiments using the finite difference method indicate 
that we can enhance the convergence rate in subcritical cases by 
extending the X-domain. The further downstream we place the 
boundary, the more the disturbance is damped near the boundary and 
this, in turn, speeds up the global convergence properties of the 
scheme. Conversely, for supercritical disturbances the further away 
the downstream boundary is placed, the larger the value of the 
disturbance amplitude becomes, and this results in a further increase 
in the number of iterations required for convergence. 

We next consider the nonlinear problem with a hump height of h = 
1.0. In this case the mean flow is still attached everywhere in the 
flowfield. The disturbance solution for this case has the same shape 
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and follows the same pattern as the previous solutions for h = 0.1 
over the entire range of Strouhal numbers considered in this study. 
However, now the amplitudes of the disturbances are increased by a 
factor of approximately 10, which indicates a linear mechanism for 
the amplification of the disturbances due to increasing the height, 
at least for this range of values of h. 

The growth rate, k^, and the wavelength, A, of the T-S waves can 
be calculated from the numerical results and compared with the values 
given by the analytical theory (Duck 1985) . The pressure p(X,t) = 

Realfp^X) e”^‘ t ) , where P X (X) is calculated numerically. Now P X (X) = 

ilcX • 

Ce + a function of X which decays algebraically as X ■ and C xs 

a constant. Therefore, (dp^/dX)/^ » ik and (d 2 p 1 /dX 2 )/p 1 « -k , 

where k = k r + ik^ and x * 2*/k r . This approximate method for 

calculating k should become more accurate when k is calculated from 

higher order derivatives of the pressure. Figures 5a, 5b show the 

values obtained at each X station for the linearized time-periodic 

spectral results for the case S Q = 1.0 calculated from the second 

derivative of the pressure. For comparison purposes, the finite 

difference results for k calculated using the first derivative 

approximation for the pressure are also shown. The two methods give 

substantially the same results for X > 10. In the region 12 < X < 

40 the values of k r are within 3% of the analytical value of k f * - 

0.522 and the values of k^ within about 10% of k^ « 0.121. The 

reason for the larger % error in k^ compared to k r is that | k^ | < 

|k r j. As is to be expected, as we move upstream of this region the 

algebraic terms become significant as the influence of the hump 
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increases and the calculated values of k , k^ show greater variation 
from their analytical values. 

Figure 6 shows the results obtained for h = 1.0, S Q = 3. The 
finite difference results were calculated using an X grid -10 s X s 
15 with AX = 0.125. The time-marching spectral method results here 
are for t * 5* (with the negative of the results plotted for 
comparison with the finite difference solution) . By this time the 
time-periodic solution is established in the region approximately X < 
13 and here there is excellent agreement with the finite difference 
results. 

Next, results for h = 5, for which the mean flow has a small 
region of separated flow, are presented. In Figure 2, the 
corresponding mean flow distributions for the pressure and wall shear 
are presented. We note that the solutions obtained by the finite 
difference method and the steady spectral method are in excellent 
agreement. There are some small oscillations in the wall shear 
results in the reversed flow region suggesting the need for a 
somewhat finer grid for the finite difference method. This effect 
was confirmed by running the finite-difference code on a somewhat 
finer grid (-10 s X s 10, 400 points) . In this case the oscillations 
noted above did indeed disappear. 

For this value of h, the finite difference method was used to 
obtain converged solutions for S Q = 1 and 1.2. Pressure and wall 
shear distributions for the unsteady flow over the X-domain, -10 s X 
s 10 is given in Figure 7 for S Q -1. 
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Converged solutions using the finite difference method for 

larger values of the Strouhal number are possible in principle; 

however, the computer time necessary to achieve them becomes 

prohibitive. Furthermore, the accuracy of the numerical solutions 

deteriorates as the streamwise step size, AX, is increased. This 

deterioration can best be seen in the wall shear results. 

Maintaining the same number of mesh points (200) and increasing AX we 

have been able to find three other "numerically converged" solutions 

for h =5 and *= 1. The wall shear distributions for these 
o 

solutions are given in Figure 8. Note that for all cases the general 
shape of the solutions do not change. However, by increasing AX the 
magnitude of the disturbance amplitude was reduced while maintaining 
its same form and location. For the largest value of AX, which 
corresponds .to an X domain of -10 s X s 90 (AX=0.50), the numerical 
breakdown is clearly seen in Figure 8c even though the general shape 
of the distribution has not changed appreciably. For this same X- 
domain, another "numerical solution" was obtained for half the step 
size by doubling the number of points to 400 and the results are 
shown in Figure 8d. Note that the numerical oscillations have 
disappeared and the same general shape of the curve remains. 

However, the magnitude of the amplitude is the same as that in Figure 
8b. There thus seems to be a trade-off between a larger X-domain to 
help on the convergence rate and a smaller step size to maintain an 
accurate solution. 

For h - 5, S ■ 1, the finite difference results calculated on 
the three grids; -10 s X s 40, AX * 0.25; -10 s X s 20, AX = 0.15; 
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and -10 s X ^ 10, AX = 0.1; giving discrepancies in the amplitude of 
the T-S wave of up to 30%, 10% and 5%, respectively, when compared 
with the time-periodic spectral method results which were calculated 
on a grid sufficiently fine and extensive to give grid independent 
results to graphical accuracy. A comparison of these results for the 
spectral method and the finite difference method using the third grid 
is also given in Figure 7. 

A comparison between the numerical computations and Goldstein's 
(1985) analytical theory for h « 1 can also be made. The 
disturbance amplitude computed is shown in Figures 9a, b for h = 1, S Q 
* 3 and h = 5, S Q = 1.2, respectively. For h = 1, the analytical 
theory and the numerical results are in good agreement over most of 
the region of interest. For smaller values of h and other values of 
' S Q , not shown here, indicate an even better agreement between the two 
approaches. Thus we can conclude that the analytical theory can be 
applied for h < 1. Conversely, for h = 5, the disturbance amplitude 
predicted by Goldstein's theory does not agree well with that 
predicted by the numerical computations. The nonlinear base flow 
results in substantial differences of the receptivity problem from 
that computed by the analytical theory. 

We next consider the physical implications of the numerical 
results, and in particular the effect of the hump height, h, on the 
disturbance wave amplitude. We choose to follow the amplitude 
variation with h of several peaks in the disturbance wall shear 
solution. These peak values are scaled by the corresponding peak 
values of the linear results: Note that the location of the peaks 


22 


does vary slightly with h. The peak amplitude dependence on h was 
studied for subcritical values of S Q , namely S Q = 1,2 and for the 
supercritical value S Q = 3. The S Q = 3 solutions were found for h s 
3.5, calculations for larger values of h were affected by a grid 
dependent rapid growth in the spectral plane solution for |k| » 1. 
This possible Rayleigh instability (see Smith 6 Bodonyi 1985, Tutty & 
Cowley 1986) is suppressed in the time-periodic spectral method and 
the finite-difference method which treats the "steady-state" 
equations to find solutions which are periodic in time. 

Consequently, solutions can be found for Is h s 5. 

In Figures 10a, b we show the results for S Q = 1 and S Q = 3, 
giving the scaled peak amplitude behavior for the wall shear r(x,0) 
with h for two different peaks; the results for S 0 = 2, not presented 
here, follow a similar pattern. 

For h < 1, the scaled disturbance amplitude depends approximately 
linearly on the hump height. Experimentally, Azin & Polyakov (1979) 
found a linear dependence of the disturbance amplitude on h, for the 
interaction of upstream propagating sound waves with thin mylar 
strips fixed on a flat plate near the lower branch of the neutral 
stability curve. For larger values of h, our results show an 
increasingly nonlinear enhancing effect on the disturbance amplitude. 
These results are confirmed by the finite difference results for S D = 
1. For example, for h * 3 the disturbance amplitude is approximately 
twice that of the linear results for the same h. The results for the 
subcritical S Q calculated using the time-periodic spectral method and 
finite difference method show that this increasingly rapid 
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enhancement of the "receptivity" continues for h ■ 4 and 5. The 
results for h = 4,5 are approximately three and six times, 
respectively, the corresponding linear results for the same h. 

In conclusion we find a linear dependence of the disturbance 
amplitude on the hump height for sufficiently small values of h. For 
moderate h(lshs3)we find an enhancement of the receptivity by 
the nonlinear effect of hump height. For large h where local flow 
separation can occur in the steady flow, we find a possible short 
wavelength instability in our time marching calculations and a 
rapidly increasing enhanced receptivity in our "steady state" 
calculations . 
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Wall shear distribution for steady flow with h ■ 5. 
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Disturbance amplitude variation versus h for 
normalized by the linearized results for h - 

S Q - 1. : minimum near X - 18.5; - - - : 

near X * 13; : linear result. 
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Figure 10b. 


Disturbance amplitude variation versus h for s n “ 3 
normalized by the linearized results for h ** l, 

S 0 - 3- : minimum near X - 9; : peak near 

x m •• : linear result. 



